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ABSTRACT 

Semi-analytic models (SAMs) are currently one of the primary tools with which to 
model statistically significant ensembles of galaxies. The underlying physical prescrip- 
tions inherent to each SAM are, in many cases, different from one another. Several 
SAMs have been applied to the dark matter merger trees extracted from the Mil- 
lennium Run, including those associated with the well-known Munich and Durham 
lineages. We compare the predicted luminosi ty dis tributions of galaxy groups using 
four publicly available SAMs (|De Lucia et al.ll2006t iBower et al]|2006t iBertone et all 
l2007tlFbnt et al. 2008), in order to explore a galactic environment in which the models 
have not been explored to the same degree as they have in the field or in rich clusters. 
We identify a characteristic "wiggle" in the group galaxy luminosity function gener- 
ated using the De Lucia et al. (2006) SAM, that is not present in the Durham-based 
models, consistent to some degree with observations. However, a comparison between 
conditional luminosity functions of groups between the models and observations of 
Yang et al. (2007) suggest that neither model is a particularly good match. The lu- 
minosity function wiggle is interpreted as the result of the two-mode AGN feedback 
implementation used in the Munich models, which itself results in flattened magnitude 
gap distribution. An associated analysis of the magnitude gap distribution between 
first- and second-ranked group galaxies shows that while the Durham models yield 
distributions with approximately equal luminosity first- and second-ranked galaxies, 
in agreement with observations, the De Lucia et al. models favours the scenario in 
which the second-ranked galaxy is approximately one magnitude fainter than the pri- 
maryespecially when the dynamic range of the mock data is limited to 3 magnitudes. 

Key words: Galaxies: luminosity function - galaxies: clusters - galaxies: formation 
- methods: N-body simulations 



1 INTRODUCTION 

In an ideal world, simulating and analysing the formation 
and evolution of galaxies within Gpc-scale cosmological vol- 
umes would be accomplished via the use of self-consistent 
gravitational N-body and hydrodynamical models. While 
such an approach is feasible in certain restricted situations, 
it remains impractical for most applications. Instead, the 
compromise most widely adopted is based upon the use of 
Semi- Analytic Models (SAMs). The merger and assembly 
histories of galaxies within SAMs are underpinned by high- 
resolution cosmological N-body simulations, at the "cost" of 



employing a posteriori "semi-analytical" treatments of the 
associated baryonic physics. 

In all the SAM models explored in this paper galaxy 
properties are derived using a range of gas infall, radiative 
cooling, re-ionisation, AGN and supernovae feedback, mor- 
phological transformation, dust and spectrophotometry pre- 
scriptions. In general, the inclusion of AGN feedback within 
the SAM reduces the luminosity and stellar mass of the 
brightest galaxies. Supernovae feedback is effective in low 
mass galaxies, where it becomes an important mechanism 
by which galactic winds are driven and star formation is 
quenched. Thus, supernova feedback leads to a reduction in 
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the number of low luminosity galaxies. Upon merging, the 
stars and cold gas of a satellite galaxy (an accreting galaxy) 
are added to the reservoir of the central galaxy (called 'Cen- 
trals', henceforth) of the parent halo. SAMs have reproduced 
a range of galactic observables, including colours, luminosi- 
ties, and mass functions. 

SAMs come in several flavours, and, although the 
codes share many similar features as outlined above, they 
also differ in the way in which certain processes, relat- 
ing to baryonic physics, are implemented (e.g., treatment 
of supernovae and AGN feedback). These lead different 
SAMs to produce different solutions to the problem of 
galaxy formation. Some of these differences have been 
explored at length in the literature, via direct compari- 
son with both empirical field and cluster galaxy luminos- 
ity fun ctions, which are t he extrema of galaxy environ- 
ment s jHatton et al.l 12003 ; iMo et al.l |2004| ; iGonzalez et al.1 



120051; iBower et al l 1200(1 ICroton et all 120061). For e xam- 
ple. iMateusl (l200Sh fo und that the IBower et alJ l|2006h and 



iDe Lucia et al. I |2006h models give different trends for the 
temporal evol ution of galaxy merger r a tes b ased on close 
pair counting. | Dfaz-Gimenez fc Mamonl (l2010h sugge st that 
the IBower et all (|2006r i and IDe Lucia et al.1 (|2006r ) mod- 
els reach different conclusions regarding the rate of chance 
alignment in low velocity dispersion compact groups. Recent 
examples of problems encountered by the SAM approach in- 
clude the excess of low -mass red galaxies, as ide ntified by 
IWeinmann et al l l|2006h . and iBaldrv et al.l (|2006h . A com- 
prehen s ive re view of the SAM approach can be found in 
iBauehl J200ri. 

What has not been explored, thus far, at least in any for- 
mal sense, is the impact of these baryonic physics prescrip- 
tions upon the resulting luminosity and stellar mass func- 
tions for the most common of environments, that of galaxy 
groups. It is to this aim that our current study is focused. 
Galaxy groups are environments where galactic evolution is 
happening at a high rate due to the low velocity dispersion 
of groups. This means means that galaxy-galaxy interactions 
are more likely than in clusters. In this paper, we examine 
the outputs of four widely-used SA Ms applied to the Millen- 
nium RunQ Springel et al.l (|2005h . in order to quantify the 
impact of baryonic physics prescriptions upon the resulting 
compact and loose group luminosity functions. Two of the 
models which we examine will be colle ctively referred to a s 
the "Durham models", b e ing th ose of IBower et al.l {2006), 
(D_B06), and iFont etaP (120081 ) (D_F08), which is an up- 
dated version of DJ306, with a more sophisticated treat- 
ment of ram pressure stripp ing. We also analyse t wo "Mu- 
nich models", being those of IDe Lucia et al.1 (|2006h (M_D06 
hereafter), and lBertone et al.1 (|2007f) (M_B07), which differs 
from M_D06 mainl y in the supernov a feedback recipes. A 
related model bv lCroton et al.l (|2006l ). of which M_D06 is a 
direct descendant, is also referred to in our study. 

After outlining our galaxy group cataloguing proce- 
dure, constructed using a classical friends-of-friends ap- 
proach (§ 2), we examine systematically the predicted dis- 



1 The simulation was carried out by the Virgo Super- 
computing Consortium at the Computing Centre of the 
Max Planck Society in Garching and recovered from 
http : //galaxy-catalogue . dur . ac . uk: 8080/Millennium/ . 



tributions of luminosity, and first-to-second rank magnitude 
gap for both compact and loose groups of galaxies, for each 
of the SAMs under consideration (§ 3). We analyse the lumi- 
nosity distribution of galaxy groups in the different models 
so that the next generation of SAMs can improve the imple- 
mentation of galaxy formation physics. 



2 MODELS 

The SAMs used in our analysis employ th e merger trees 
assoc iated with the Millennium Simulation (|Springel et al.l 
2005) a large N-body simulation corresponding to a sig- 
nificant volume of the visible Univer se, and generated u s- 
ing the WMAP Year 1 cosmology <|Spergel et al.l l2003h FI 
The simulation used 2160 3 particles in a periodic box of 
side length 500ft" 1 Mpc, gravitational softening of 5h _1 kpc, 
and individual particle masses of 8.6x10 s M ; 64 outputs 
exist within the Millennium database, ranging from red- 
shift 2=127 to z—0. The simulation was p ost-processed us- 
ing a Friends-of-Friends (FoF) algorithm dGeller fc Huchral 
1983) , in order to ide ntify density peaks, and then SUBFIND 
1 Springel et al.1 2001) was employed to identify substructure 
and split spuriously joined haloes. This information was then 
used to build merger trees for the dark matter haloes, onto 
which the SAMs are "mapped" . 

Before embarking on a discussion of the predicted lu- 
minosity functions resulting from the use of the aforemen- 
tioned SAMs applied to the Millennium merger trees, it is 
important to summarise briefly the defining characteristics 
associated with each of the primary SAMs employed here. 
We will highlight the different ways in which the codes create 
merger trees, the way in which galaxy positions are defined, 
the implementation of satellite disruption and accretion, and 
the way in which supernova and AGN feedback are imple- 
mented. 



2.1 Durham models (|Bower et al.ll2006l : 
IFont et al.ll2008h 

In the Durham models, m erger trees are pro duced in a man- 



ner which follows that of iHellv et al.l 1 2003h. and the prop - 
erties of these trees are described in Harker et al. I (|2006h . 
These models account for ostensibly separate haloes which 
are joined by a bridge of dark matter and hence can be erro- 
neously put in a single halo by FoF algorithms, and also ac- 
count for haloes which are only temporarily joined. Account- 
ing for these effects results in a halo catalogue containing 
more haloes than in the original FoF catalogue. The merger 
trees are then constructed from these catalogues by follow- 
ing subhaloes from early times to late times. We note that 
the merger tre e s wer e constructed independently of those in 
ISpringel et all (|2005l ). 

The merging of galaxies, and lifetime of s atellite galax- 
ies, ar e derived using the method presented in lBenson et al.l 
(2002), which i s considerably m ore sophisticated than the 
method used in ICole et al.l l|2000ft . When dark haloes merge, 
a new combined dark halo is formed, and the largest of 
the galaxies contained within is assumed to be the central 
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galaxy, whilst all other galaxies within the halo are satel- 
lites. Satellites are then evolved under the combined effects 
of dynamical friction, and tidal stripping. These effects are 
modelled analytically. The initial orbital energy, and angular 
momentum, of the satellite upon merging are specified. The 
orbital energy is set using a constant value of r c (E) / R V i r 
= 0.5, representative of the median binding energy of satel- 
lites, while the orbital ellipticity, is chosen to be between 0.1 
and 1.0 at random. Given these parameters, the apocentric 
distance is found, and the orbit equations are integrated at 
that point. The host and satellite haloes are all assumed to 
have NFW profiles, while galaxies are modelled as a disc 
plus spheroid. The satellite galaxies plus halo are then ad- 
vanced by calculating the combined gravitation forces of the 
host and satellite haloes, as well as the effects of dynami- 
cal friction, calculated using Chandrasekhar's formula. The 
code keeps track of tidal stripping, to remove mass from 
the satellites. The new halo mass is then used for the next 
iteration of the orbit equations. 

This integration of the orbital equations continues until 
one of three conditions are met: (i) the final redshift, (ii) 
the host merges with a new larger halo, in which case the 
satellite becomes a part of the new halo and has new orbital 
parameters assigned, or (iii) the satellite merges with the 
central galaxy. In this last case, merging takes place when 
the orbital radius falls below R m erge which is the sum of 
the half mass radius of the central and satellite galaxies. 
The mass of the merged satellite is then added to the cen- 
tr al galaxy. The m erging times match largely match those 
of lCole et ai] f2000). however some satellites have very long 
merging times, as they loose a great deal of mass through 
tidal stripping, meaning that dynamical friction forces be- 
come very weak. 

The Durham models relate supernova reheating directly 
to the cir c ular velocity of the galaxy disk according to 
ICole et all (|200Ch : 



M reheat oc V'M* 



(1) 



where M re heat is the rate of change of mass of reheated gas, 
Vdisk is the disk circular velocity, M, is the time deriva- 
tive of stellar mass and M e j ect is the change in the mass of 
ejected gas. In haloes with a shallow potential, this has the 
effect of reducing the amount of cold gas available to form 
stars, by heating the gas back into the hot gas reservoir. 
The hot gas is dominated by ejection for low mass haloes, 
and by reheating without ejection for large haloes. For low 
mass galaxies, supernova feedback is an important mecha- 
nism by which galactic winds are driven and star formation 
is quenched. 

The Durham models implement AGN feedback is such a 
way as to regulate the cooling of hot gas. In large haloes with 
large Eddington Luminosity, the AGN feedback is assumed 
to balance heating and cooling, thus truncating star forma- 
tion. This prevents the formation of over-luminous galax- 
ies. While feedback is active, the black hole is assumed to 
grow proportionally to the cooling luminosity, and gas is 
accreted due to disk instabilities. The model assumes quasi- 
hydrostatic cooling for AGN active galaxies, and has a strict 
transition between AGN "active" and "inactive" phases. 
AGN feedback becomes efficient in galaxies of mass greater 
than ~2xl0 11 /i" 1 M . 

The essential difference in the D_B06 and D_F08 models 



is the implementation of ram-pressure stripping of the hot 
gas. In the D_B06 model, along with both Munich models, 
the hot gas is instantaneously stripped when it enters a halo 
already containing a central galaxy. In the D_F08 model this 
process happens gradually and depends on the orbit of the 
galaxy. This has the effect of reducing the population of faint 
red galaxies. 



2.2 Munich models dDe Lucia et al.| [2006; 
iBertone et al.ll2007h 



The Munich merger trees (|Springel et al l l2005h used in 
M_D06 and MJ307 follow the positions of subhaloes for as 
long as they can be identified. Identification is possible only 
when the number of particles bound to a subhalo exceeds the 
minimum number of particles set by SUBFIND. The trees 
are constructed by following the most bound halo particles 
and searching for the descendant halo in the next output. 

One of the key differences between the Munich mod- 
els and those of Durham is that the Munich models explic- 
itly follow dark matter haloes even after they are accreted 
onto larger systems, allowing the dynamics of satellite galax- 
ies residing in the infalling haloes to be followed until the 
dark matter substructure is destroyed. The galaxy position 
is calculated by assigning the galaxy to the most bound par- 
ticle of a (sub)halo at each time step. This is done until the 
(sub)halo is no longer identifiable, whereupon the galaxy 
is assigned to the most bound particle of the (sub) halo at 
the last time the (sub)halo could be identified. An analytic 
countdown to galaxy merging begins when the satellite sub- 
halo can no longer be identified, and resets if the parent 
halo undergoes a major merger. Thus, in the Munich mod- 
els, the lifetime of galaxies in groups depends on the amount 
of time the (sub)halo finder can identify the subhalo, plus 
th e analytic countdow n. The analytic merging follows that 
of ICroton et"afl j200rj ): 



1.17 



V H r 2 aat 



Gm 3at In (l + ^)' 



(2) 



where hi (l + M - - j is the coulomb logarithm, Vh is the halo 
circular velocity, r sat is the distance of the subhalo from the 
halo centre at the time it is last identified, m sa t is the mass 
of the satellite dark halo at the time it was last identified 
and Mh is the halo mass. 

In M_D06, the amount of reheated (by supernovae) cold 
gas is proportional to the stellar mass, and the mass ejected 
from the halo is inversel y proportional to t he host halo's 
circular velocity squared (|Croton et al.|[200r? ); 



M reheat OC V vir M* 



(3) 



The important parameters, Vvir and M», have the same re- 
lationship in the Durham models. In the Munich models, 
the reheated mass is proportional to the mass of the halo, 
while in the Durham model it is proportional to the disk 
mass. M_B07 adopts a more sophisticated treatment of su- 
pernova feedback. Rather than simply parametrising the ef- 
fect of supernovae feedback, the M_B07 model follows the 
dynamical evolution of the wind as an adiabatic expansion 
followed by snowploughing. This implementation has the ef- 
fect of increasing the luminosity of the brightest galaxies. As 
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also occurs with the Durham models, some of the gas will 
be ejected from low mass haloes in both Munich models. 

In the Munich models, a two-mode formalism is adopted 
for AGN, wherein a high-energy, or "quasar" mode oc- 
curs subsequent to mergers, and a constant low-energy "ra- 
dio" mode suppresses cooling flows due to t he interaction 
betwe en the gas and the central black hole ICroton et al.l 
(2006). In the quasar model, accretion of gas onto the black 
hole peaks at z ~3, while the radio mode reaches a plateau 
at z ~2. AGN feedback is assumed to be efficient only in 
massive haloes, with supernova feedback being more domi- 
nant in lower-mass haloes. 

The properties of groups in the Munich models, MJD06 
and M_B07, are similar to one another in many cases, as 
are the properties of the two Durham models, D_B06 and 
D_F08. Thus, in some analysis, we just discuss the MJD06 
and DJ306 models, as representatives of their model's lin- 
eage. We only discuss the results based on their descendants, 
M_B07 and D_F08, when they show significantly different 
behaviour from M_D06 and D_B06. 



3 GROUPS 

In order to construct a statistically significant (and represen- 
tative) galaxy group catalogue, we have worked with a set 
of sub-samples of the Millennium Simulation, amounting to 
~3% of the available volume, - specifically, 64 boxes of side 
length 125/i _1 Mpc drawn from the database. Our results 
are robust to the arbitrary selection of the box, having been 
tested a posteriori on alternate boxes of equal size. A lumi- 
nosity limit of M r =— 17 in the SDSS r-band was imposed. 
At lower luminosity the effect of the limited mass resolu- 
tion of the N-body background effects the completeness of 
the sample. We identify galaxy groups as o verdensties in the 
galax y population using an FoF algorithm (|Geller fc Huchral 
1 19831 ). No maximum number of members is set but we re- 
quire that at least four galaxies are linked in order to define 
a group. Although t his removes gro ups such as the Local 
group, it follows the iHicksonl (|l982l ) definition of compact 
groups. 

We first construct a "loose group" (LG) catalogue using 
a linking length of 0.2 times the mean inter-particle sepa- 
ration, (or in this case, inter-galactic separation). We made 
this choice by assuming that the galaxies follow the dark 
matter. This corresponds to a co-moving linking length of 
^OOh' 1 kpc. To examine the effects of density we also 
define two "compact group" catalogues, - Compact (CG), 
and Very Compact (vCG), Groups - using co- moving link- 
ing lengths of 150/i _1 kpc and 50/i _1 kpc, respectively. The 
C G linking length of \ 5 0h~ 1 kpc is similar to that advised 
bv lMcConnachie et al.l j200§), based upon their 3D linking 
length analysis from mock catalogues of Hickson compact 
groups based on M_D06. The vCG linking len gth is compa- 
rable to the projected linkin g lengt h used bv lBarton et al.l 
1 19961 ) and lAllam fc Tucked (|200fj| ) to identify groups, a 
value arrived at by calibrating to the iHickson et all (|l992f ) 
catalogue using projected galaxy separations. 

We note that the CG galaxies are, by necessity of the 
group finding algorithm, subsets of the LG catalogue, in that 
every galaxy assembled into a group at short linking length, 
must be part of a group with a larger linking length. Our 



Table 1. The number density of field galaxies, and LG, CG and 
vCG for the D_B06,M_D06, M.B07 and D.F08 models . 





D.B06 


D.F08 


M.D06 


M.B07 




ft 3 Mpc~ 3 


/i 3 Mpc -3 


/i 3 Mpc -3 


/i 3 Mpc -3 


Field 


6.5 x 10 2 


6.6 x 10 2 


6.3 x 10 2 


5.0 x 10 2 


LG 


2.1 x 10 3 


2.2 x 10 3 


2.0 x 10 3 


1.6 x 10 3 


CG 


1.3 x 10 3 


1.4 x 10 3 


8.9 x 10 4 


8.7 x 10 4 


vCG 


5.4 x 10 4 


6.6 x 10 4 


1.2 x 10 4 


1.4 x 10 4 



catalogues also contains clusters and cluster cores, a point 
to which we return shortly. The physical interpretation of 
the linking length variation and its impact upon resulting 
galaxy distribution is non-trivial. The FoF algorithm essen- 
tially probes deeper into the potential well at shorter linking 
lengths, selecting only galaxies closer to the cluster/group 
core. These galaxies are generally old, and have sunk deeper 
into the cluster potential, or they are galaxies near their 
respective orbital peri-centre. 

The algorithm does not only extract the inner region of 
groups. Galaxies in the outskirts of clusters that happen to 
be close to one another can also be linked together. This may 
either be due to a temporary alignment of galaxies that are 
simply 'passing through' that region, or because the galax- 
ies were in close proximity before they entered the cluster, 
and have not yet been disbursed by tidal forces. These pe- 
ripheral groups are essentially cluster substructure, and are 
a natural part of the analysis. Peripheral groups represent a 
small fraction of LGs, but represent ~20% of CGs and vCGs. 
Also worth noting, the linking length adopted for LGs can 
include galaxies in haloes beyond the limit of the dark mat- 
ter halo which the majority of the group members occupy. 
This means that our LG catalogue is not exactly equivalent 
to groups that are defined by the halo occupation of galaxies. 

3.1 Merging Timescales 

We compare the merging histories and timescales of haloes 
and galaxies within the models. Any differences will be im- 
portant in determining the origin of various group proper- 
ties. To compare the lifetime of satellites we look at a) the 
lifetimes of satellites that have merged, and thus contribute 
to the mass and Luminosity of the Central Galaxy, and b) 
the lifetime of satellites that have not yet merged, and are 
hence satellites at z = 0, and contribute to the Group cata- 
logues. 

Figure [T] shows the distribution of merging times, de- 
fined as the time between a satellite first entering a halo 
and being totally merged with the central galaxy. We se- 
lect the same 19 millennium simulation clusters from each 
model to make a fair comparison between the codes and do 
not put any limit on the luminosity of infalling satellites. 
The distribution of merging times shows little difference be- 
tween the models. The average M_D06 satellite has lasted 
5.8 Gyr, while D_B06 satellites last 6.2 Gyr, with standard 
deviations of ~ 1. We note, however, that twice as many 
satellites have merged in the Munich model (317) than in 
the Durham model (152). 

This result seems a little contradictory, as similar merg- 
ing timescales should result in similar numbers of mergers. 
We look to the number of satellites which do not merge in 
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Figure 1. Shows the distribution of galaxy merging times. The 
y-axis is number of galaxies in each bin normalised by the total 
number of galaxy mergers. The black line is the D_B06 model and 
the red dotted line is for the MJ306 model. Errors are Poisson. 
The number of merging galaxies in each model is quoted in the 
panel. 




Figure 2. Panel A: Shows the number of galaxies which merge 
with the halo central galaxy during the simulation against the 
time the galaxy entered the host halo. Panel B: Shows the infall 
time of those galaxies which are still satellites at z=0. Panel C: 
is the sum of the plots in Panel A and Panel B. The solid black 
line shows the D_B06 result and the red dashed line is for M_D06. 



order to reconcile this. Figure [2] shows the time of infall for 
the satellites into the host halo, for satellites that merge 
(Panel A), satellites that do not merge but rather remain 
as satellites at z=0 (Panel B), and all satellites (Panel C). 
Again, the fact that more satellites merge in M_D06 com- 
pared with DJ306, can be seen in panel A, but it can be 
seen that the difference is dominated by satellites which fall 
in to the host halo at early times. By contrast, the num- 
ber of satellites which do not merge at all is significantly 
larger for DJ306, compared with M_D06. The difference is 
greatest for satellites which accrete early. This significant 
population of satellites which do not merge by z=0 can be 
expected to affect the properties of groups which we present 
in the remainder of this study, both because of the lower 
number of satellites "feeding" the central galaxy in DJ306, 
and also the larger numbers of satellites which survive to be 
included in the group catalogues. Panel C shows the infall 
time of all satellites from our sample of 19 clusters, regard- 
less of whether they merged, and shows that the differences 
between the models being relatively minor. 

There are more haloes falling into the Durham models 
and fewer galaxy mergers. The galaxies which merge have 
the same lifetime in each model. If we look at the time of 
infall of galaxies which merge, and those that remain in 
the cluster at z~0, (Fig. [2}, there is a substantial popu- 
lation of galaxies which never merge. In the Durham mod- 
els more galaxies survive than in the Munich models (panel 
(B) of Fig. 2). Further, we also see that more (of the early 
accreted) galaxies merge with the central object in the Mu- 
nich model (panel (A) of Fig.2).At the earliest times in the 
M_D06 model almost all the galaxies have merged, while in 
the DJ306 model they have not. Thus, the 'maximum' life- 
time of satellite galaxies is lOGyr in M_D06 but > 13Gyr in 
D_B06. 



While not shown, our examination of the merger trees 
also found that in the DJ306 model a greater number of 
haloes merge with the 19 clusters. This is because subhaloes 
are not followed when constructing the merger trees in the 
Durham models. In the Munich models there is a delay in 
the halo merger time relative to the Durham models became 
infalling haloes are able to enter a subhalo stage before they 
are considered merged. When the (sub) halo in each model 
is no longer identifiable, the SAMs are used to calculate the 
lifetimes of galaxies. In the Durham models there appears 
to be an extremally long lived population of satellites which 
are not present in the Munich 



4 OBSERVATIONS 

We use t h e ob se rvational lo ose group c atalogues of 

iTaeo et all l|200Sl), lYaiig et all (I200S ). and iTucker et all 
(|200fj|) . and the lAllam fc Tuckerl (|2000h compact group cat- 
alogue. These group catalogues take large redshift surveys 
of gala xies and use a FoF algorithm to assemble them into 
groups. [Tago et all (|2008l . SDSS5.T08) implemented a stan- 
dard FoF algorithm which scales according to the distance, 
and a pplied it to the SDSS DR5 l|Adelman-McCarthv et al I 
l2007h . The initial linking length is 0. 25/^^ 00 in projec- 
tion, and 250fcms _1 along the line of sight. I Tucker et al.l 
l|200d. LCRS TOO ) catalo gue uses the Las Campanas Red- 
shift Survey (jShectman et alJll996t ). and applies a fiducial 
linking length of 0.715/i _1 Mpc and 500 km s~ 1 , scaled from 
this value at z = 0.1 a nd rising to 1.8 times this at z = 1.7. 
lAllam fc Tuckerl (|2000l . LCRS_A00) uses the same catalogue 
as LCRS.TOO but a sh orter linking leng th of 0.05/i _1 Mpc 
and 500/cms" 1 . Finally. lYaiig et al l (|2008l . SDSS4.Y08) uses 
a more complicated iterative approa ch, which, nevert heless, 
includes a FoF algorithm at its core (|Yang et al.ll2005T ). This 
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algori thm applied to SDSS DR4 (|Adelman-McCarthv et all 
l2006h . 

Our synthetic group catalogues are compared directly 
with these empirical datasets wherever such comparisons are 
possible and meaningful within our analysis, bearing in mind 
limitations in the models, and in the observations. Where 
appropriate, "cuts" are made in the synthetic catalogues to 
mimic observations. 



5 RESULTS 

5.1 Density of Groups in the Models 

The difference in the proportion of galaxies in groups should 
provide a key diagnostic for comparing and discriminating 
between the four SAMs, relative to empirical data. Within 
the Millennium Simulation volume at z~0, the four SAMs 
under consideration yield differing numbers of galaxies, de- 
spite being built upon the same underlying dark matter dis- 
tribution. It is useful to review the respective galaxy num- 
bers and distributions for these models. Restricting the dis- 
cussion to the relevant sampling criteria (i.e., M r < — 17) the 
M_D06 and M_B07, and D_B06 and D_F08 models have a 
galaxy number density given in Table [T] The number den- 
sity of galaxies in the field for each model is quite similar, 
as well as in LGs. However, the Durham models have sig- 
nificantly denser galaxy populations within CGs and vCGs. 
The physics which the different models employ to account 
for group environments appear to have significantly altered 
the nature of compact groups. 

The relative proportions galaxies that are classified as 
being members of groups, along with the average group rich- 
ness, are listed in Table I5TT1 In this instance we define group 
richness as simply being the number of galaxies in a group. 

The percentage of galaxies associated with LGs (and 
the numbers of galaxies per loose group) is comparable 
between three of the four SAM variants, with the M_D06 
model showing approximately 6 percent fewer groups than 
the DJ306 model. The models diverge increasingly with de- 
creasing linking length, with the Munich models, M_D06 and 
M_B07, having 5 to 6 times fewer vCGs than in the Durham 
models. The M_D06 model produces noticeably fewer rich 
groups at all linking lengths compared with the other Mu- 
nich model, M_B07. This indicates that the implementation 
of supernovae has a significant effect on the richness of group 
catalogues. The Durham models, DJ306 and D_F08, have 
slightly richer loose groups than the Munich models. This is 
associated with the smaller populations of medium bright- 
ness red galaxies in the Munich models, which may be the 
result of the differences in the creation of halo catalogues, 
the tracing of subhalo mergers, or due to radio mode AGN 
feed back. 

iMcConnachie et al.l |2008l . 120091 ) compared compact 
groups in mock redshift catalogues to SDSS DR6 observa- 
tions, and concluded that the M_D06 SAM overproduces 
CGs by ~50%. By extension, Table 15.11 indicates that the 
DJ306 and D_F08 models result in an even more dramatic 
"overproduction" of CGs (by an order-of- magnitude). Thus, 
in this regime, none of the models are good fits to the em- 
pirical data, and the Durham models are particularly poor. 

Figure [3] shows the galaxy distribution of galaxies with 



Table 2. Top: Percentage of galaxies in FoF groups for LG, CG 
and vCG in each of the four SAMs; Bottom: Mean group richness 
for groups in the four models 





LG 


CG 


vCG 


Bower (D_B06) 


44 


18 


5.5 


Font (D.F08) 


41 


19 


6.0 


De Lucia (M.D06) 


38 


10 


0.9 


Bertone (M.B07) 


43 


14 


1.3 


Bower (D.B06) 


13.5 


9.4 


6.6 


Font (D.F08) 


13.1 
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Dc Lucia (M.D06) 
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Bertone (M.B07) 
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Figure 3. Left Panel: The number of galaxies in LGs against 
distance from the centre, normalised by the total number of galax- 
ies. The solid line is for D_B06 LGs and the dashed line is for 
M_D06 LGs. Right Panel: The normalised number of galaxies 
against the radius scaled by the virial radius of the host.The red 
line shows SDSS4_Y08 LGs for the 60 most massive clusters. Er- 
rors are Poisson. There are 2294, 2022 and 4688 in the D.B06, 
M.D06 and SDSS4.Y08 datasets respectively. 
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Figure 4. The distribution of nearest neighbour pair separations 
for galaxies in LGs for the three models, the solid line is D_B06, 
the dashed line is M_D06, and the dotted line is DJ307. The plot 
shows that D_B06 galaxies tend to be closer together on average. 
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Figure 5. Derived SDSS r-band luminosity functions for galaxy 
groups constructed with the SAMs described herein: "All" refers 
to the global galaxy luminosity function; "LG" to loose groups; 
"CG" to compact groups; "vCG" to very compact groups. The 
solid/dashed/dotted lines refer to the D_B06/M_D06/M_B07 
models, respectively. 



radius for the LGs of M_D06 and D_B06 models. The left 
panel shows the number of satellites versus radius in Mpc. 
The right panel of Fig.[3]shows the same plot with the radius 
scaled by the virial radius. Both panels show that there are 
more satellites in the inner regions of the DJ306 LGs, rela- 
tive to the outer regions, compared with the M_D06 model. 
The red line in the right panel shows the 60 most massive 
SDSS4_Y08 observed Loose Groups, making them similar in 
mass range to the population used in making the simulation 
plots. This indicates that both models have a radial distribu- 
tion of satellites which is more centrally concentrated than 
the observations, with the problem being particularly acute 
in the D_B06 model. 

This more concentrated distribution of satellites is also 
reflected in the distribution of nearest neighbour pair sepa- 
ration distances for LG members, seen in Fig. [4] It can be 
seen that the D J306 model the nearest neighbour separation 
is smaller than for the M_D06 and M_B07. We note that the 
D_F08 model groups have very similar properties to DJ306 
and is thus not included in these figures. 



5.2 Luminosity Functions 

We now plot the group luminosity functions. Figure [5] shows 
how the luminosity function changes in the different group 
environments for the SAM models. It is unsurprising that 
the global luminosity functions of all the models are similar, 
because the semi-analytic models are design ed to replicate 
the sa me observational luminosity function of Bl anton et al.l 
(2003). At smaller linking lengths (moving from top to bot- 
tom in the figure) there is a decreasing number of galaxies in 
all SAMS, a fact which is more dramatic in the Munich vari- 
ants, as noted in § 15.11 This provides further evidence that 



the Durham model galaxies are more centrally concentrated 
than those of the Munich variants. 

The second feature of Fig. [S] is the relative dearth 
of intermediate luminosity (— 21<M r < — 18) galaxies in the 
M_D06 catalogues (in relation to a simple Schechter (1976) 
function). This is manifest in the "wiggle" or "dip" seen in 
the M_D06 group luminosity function. This feature is not 
present in the other models. This wiggle becomes more ap- 
parent at shorter linking lengths, (i.e., CGs and vCGs). The 
M_B07 SAM, which employs the same AGN feedback pre- 
scription as M_D06, shows no such feature in the luminosity 
function. Weinmann et al. (2006; fig 3) show a similar wiggle 
in the luminosity function of groups in particular mass bins, 
i.e. the conditional luminosity function, using the Croton 
et al. (2006) SAM, which is a close "cousin" to the M_D06 
SAM employed here. 

In order to better understand the origin of the shape 
of these luminosity functions, we decompose the luminos- 
ity functions of the three primary SAMs (D_B06, M_D06, 
M_B07) and the three primary linking lengths (LG, CG, 
vCG) under consideration here, into central galaxies and 
satellites (Fig. [6}. The global galaxy luminosity function 
is also shown for reference. The M_D06 (middle row), and 
M_B07 (bottom row) centrals show a Gaussian distribution, 
while the DJ306 model, (top row), centrals are better de- 
scribed by a Schechter function. There is an increase in the 
number of faint central galaxies in LGs compared with CGs 
and vCGs, which alters the shape of the central luminosity 
function as we move to larger linking lengths. All satellites 
are distributed roughly according to Schechter functions. 

The wiggle in the M_D06 luminosity function is due to 
the combined effect of (i) a general lack of satellites, and (ii) 
the "peaked" nature of the central galaxy luminosity distri- 
bution. Relative to the DJ306 SAM prediction, this distri- 
bution of centrals is very narrow, without the low-luminosity 
tail associated with the DJ306 model. This effect is partic- 
ularly apparent at the shortest linking lengths, i.e. in vCGs. 
The greater total number of vCGs in the Durham models 
compared to the Munich models is also clear in the vCG 
panels, again reflecting their cental concentration. 

The contrast between the M_D06 and M_B07 models 
is of particular interest because they use the same AGN 
feedback implementation, but differ in their choice of super- 
novae feedback. The more sophisticated, and more effective 
SNe model of M_B07 makes the luminosity distribution of 
centrals wider by making a tail towards fainter magnitudes 
that is not present in M_D06. The most luminous galax- 
ies are also more luminous in MJ307 than in M_D06. The 
M_B07 treatment leads to a more significant population of 
intermediate luminosity satellites with a shallower luminos- 
ity function, thus smoothing the wiggle. A lower number 
of low-luminosity satellite galaxies in such a shallow lumi- 
nosity function is acc ompanied by "over-l uminous" massive 
(luminous) galaxies (|Bertone et al.ll2007l ). M_B07 shows a 
decrease in the number of low luminosity satellites and an 
increase in the number of intermediate luminosity galax- 
ies, reflecting a SN feedback in M_B07 which is stronger in 
dwarfs and weaker in large haloes. The increase in high lu- 
minosity g alaxies and decrease in low luminosity galaxies 
is noted in iBertone et al] (|2007T ) which they suggest could 
be solved by increasing the supernova feedback time or in- 
creasing the effect of AGN feedback. It also produces a wider 



8 0. N. Snaith et al. 



-5 



D_B06 (LG) 


D_B06 (CG) 


D_B06 (vCG) 


r 


■ r """""" 


■ /C^- — 

r - 


M_D06 (LG) 


M_D06 (CG) ^ 
/ 


M_D06 (vCG) 






/ y^P^ ' 


■ // 

-jLui 1 1 1 1 1 : 


ft 




M_B07 (LG) 


M_B07 (CG) 


M_B07 (vCG) 
















//•' / 
//' / 






A ' 






yr / 

-■- -■- ■'- 





-23-22-21-20-19-18 -23-22-21-20-19-18 -23-22-21-20-19-18 

sdss_r 



Figure 6. Group luminosity functions for the three primary SAMs under consideration (D_B06, M_D06, M_B07). In each panel, the 
dotted lines correspond to central galaxies and the dashed lines are for the satellites; the thick solid lines represent the group luminosity 
function and the thin solid line defines the global galaxy luminosity function for reference. The legend D_B06, M_D06, M_B07 refers to 
the model as defined above and the number in brackets refers to the linking length in h — x kpc, so the first column is for LGs, the centre 
column is for CGs and the right column is for vCGs. 



central galaxy distribution by not only increasing the num- 
ber of very bright galaxies but also the population of dim 
centrals. 



5.3 Brightest Group Galaxy 

We define three types of identified groups. The first type, 
bright central groups, are those where the brightest galaxy 
is also the central. The second type, peripheral groups, are 
those without a central galaxy and, therefore, the bright- 
est galaxy is a satellite. The third type of group, dim cen- 
trals, are those where the central galaxy is not the brightest. 
Except in the D_F08 model, the central galaxy is the only 
one with hot gas and it acquires all the hot gas from in- 
falling satellites. The central galaxy is also the only galaxy 
that experiences mergers and grows hierarchically. Table 
3 shows the populations of these three types for LG, CG 
and vCG, the format being, groups with bright central 
groups/peripheral groups/dim central groups. 

There is a difference between the three main models for 
LGs. The Durham models have a smaller fraction of Bright 
Central Groups relative to dim centrals, compared with the 
Munich models. This discrepancy continues to smaller link- 
ing lengths. This likely driven by the greater number of 
mergers in the Munich models, feeding the growth of the 
central galaxy. The more sophisticated supernova feedback 



Table 3. Proportion of galaxy group types where the first number 
is the proportion of bright central, peripheral and dim central 
groups. 



Model 


LG 


CG 


vCG 


D.B06 


79/00/21 


63/12/25 


64/15/21 


M.D06 


90/01/09 


74/21/05 


75/20/05 


M.B07 


82/02/16 


70/22/08 


71/23/06 



of the M_B07 model has decreased the fraction of BGGs. 
The Munich models also have more peripheral groups at all 
linking lengths 

In Fig. [7] the luminosity functions of the first ranked 
galaxies (i.e. the brightest galaxy in each group) of our 
groups have been plotted. These have then been decomposed 
by group type, with distribution of first ranked galaxies in 
bright central groups, the first ranked galaxies in peripheral 
groups and the first ranked galaxies in dim central groups. 
It can be seen that, in the Munich SAMs, for the denser 
groups, there is a large difference between the shape of the 
LF of the brightest galaxies in central and peripheral groups. 
The difference is most extreme for the M_D06 model, where 
the low magnitude tail is due, almost entirely, to peripheral 
groups. On the contrary, in the DJ306 model the distribu- 
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Figure 7. Luminosity function of all first ranked galaxies (the brightest galaxy in each group, thin solid line); distribution of first ranked 
galaxies in bright central groups (thick solid line); the first ranked galaxies in peripheral groups (dotted line) and the first ranked galaxies 
in dim central groups (dashed line). 



tion of groups is not particularly different for the different 
group types. 



5.4 Halo-based Groups 

In this section groups consist of galaxies which lie within the 
same dark matter halo, whose extent is determined by using 
a density contour defined bya dark matter particle separa- 
tion of 0.2 times the mean inter-particleseparation. This is 
not to be confused with the linking length used todefined 
groups which acts on galaxies rather than dark matter par- 
ticles. Groups determined in this manner differ from those 
determined by using FoFalgorithms. Even for LGs there is 
not a one-to-one correspondence between halo groups and 
FoF groups. The limit on the minimum number of galaxies 
used to define a group remains at four. 

We examine the conditional luminosity functions in 
three different group mass bins. These conditional luminos- 
ity functions are plotted in Fig. for the D_B06, M_D06, 
and M_B07 SAMs. The luminosity functions are further 
separated into centrals (middle row) and satellites (bottom 
row). 

The 10 13_14 M© and 1O 14 - 15 M mass bins show little 
evidence for the "wiggle" (for any of the SAMs) . The "wig- 
gle" in the M_D06 model, and to a lesser extent the MJ307 
model, is particularly prominent in the lowest mass bin (top 
right panel), where a number of physical processes become 
relevant. Bower et al. (2006) points out that at ~ 2x 10 11 M Q 



the cooling rate exceeds the free-fall rate and the halo is no 
longer in hydrostatic equilibrium. This has repercuss ions for 
the e ffectiveness of feedback from the central source dBinnevI 
|2004 ) and is used in the Durham paper to explain the break 
in the luminosity function. Our results indicate that in the 
Munich models the processes occurring in the mass range 
may have other repercussions. 

For the highest mass bin we can see that the satellite 
luminosity function is steepest for the M_D06 SAM and shal- 
lowest for the DJ306 model. The characteristic luminosity 
at the "knee" of the Schechter function (M») is lowest for 
the DJ306 model. However, this bin has only a small effect 
on the "wiggle" alluded to in Section 5.2, which occurs at 
lower luminosity than the wiggle seen in this particular mass 
bin. In the 10 12_13 Mq bin, where the wiggle is prominent 
in the M_D06 model, the satellite distribution is fairly steep 
and the central galaxy luminosity function relatively nar- 
row and bright. By contrast the D_B06 model has a much 
broader central galaxy luminosity function, indicating a ten- 
dency to produce significantly more low-luminosity centrals 
compared to M_D06. Again, the unmerged satellites in the 
Durham models mean less feeding of the central galaxy re- 
sulting in this tail to lower luminosities. 

In order to compare this with observations we over plot 
the results of M_D06 and D J306 onto the conditional lumi- 
nosity functions provided by SDSS4_Y08. IWeinmann et all 
( 2006) note that the method used by SDSS4_Y08, presented 
in Yang et al. (2005), artificially narrows the central galaxy 
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Figure 8. The conditional luminosity functions in three mass bins which correspond roughly to clusters, large groups and small groups, 
for halo-based groups for the three primary SAMs employed here, at redshift 2=0. The top row shows the distributions of all halo 
members, the second row show central galaxies and the bottom row shows satellites. The solid, dashed and dotted lines are for the 
D.B06, M.D06 and M.B07 models respectively. 




Figure 9. The conditional luminosity function of the models, D_B06 (solid black line) and M_D06 (dashed line) plotted over the CLFs 
of SDSS4_Y08 (shown in red). The number quoted being the centre of each mass bin and the box width being ~ 0.31 mag. In panel 'a' 
the top pair of lines are the 10 14 ' 58 Af Q cut, the next is the 10 13 ' 94 M Q cut etc, down to 10 12 16 M 
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luminosity function. This is because their iterative technique 
uses the brightest galaxy luminosity in the derivation of the 
group halo mass, while in the models there is no such direct 
linking of mass and luminosity. However, the difference is 
not enough to affect our comparison. SDSS4_Y08 shows the 
CLFs for groups using SDSS DR4 galaxies, and are fit with 
modified Schechter and Gaussian functions to the satellite 
and central galaxy luminosity distributions respectively. The 
functional forms of the fits are, 



$ cen (L|M) = 



exp 



(logL-logL c ) 2 



2ai 



$sat{L\M) = 



exp 



(*)' 



(4) 



(5) 



where, L is the luminosity, L c is the mean position of the 
Gaussian, a c is the width of the Gaussian, (j>* is the nor- 
malisation of the modified Schechter function, is the low 
mass slope and log(L*) = log(L*) — 0.25 and is the position 
of the knee of the modified Schechter function. 

SDSS4_Y08 provide the best fit parameters, which we 
compare to the M_D06 and DJ306 models using the same 



mass bins as SDSS4_Y08 in Fig. |9(a)| and Fig. |9(b)| where 
Panel 'a' shows the satellite galaxy distribution and Panel 
'b' the central galaxy distribution. The highest mass bin in 
the models corresponds most closely with the observations. 

The shape of the conditional luminosity functions for 
the satellites are considerably different for the CLFs. At all 
masses, there are far fewer low luminosity satellites in the 
models than than in the observations. The discrepancy is 
less severe for the M_D06 model for high mass clusters. The 
difference between observations and models for the satellite 
luminosity function is larger as we go to lower mass. As we 
go to 10 12 ' 75 Mq and below, it is the DJ307 groups which 
have more low luminosity satellite galaxies. 

The central galaxy conditional luminosity functions also 
differ significantly between models and observations. The 
M_D06 model shows that the low mass (10 1216 Mq) group 
centrals peak in the same place as the observations, while 
the mass bins are somewhat displaced, both toward lower 
luminosities (10 12 ' 75 M and 10 13 ' 3 M Q ) and higher lumi- 
nosities (10 13 ' 94 M and 10 14 ' 58 Mq). The M_D06 model is 
also broader than the observations. These discrepancies are 
even greater for the DJ306 groups, which even wider than 
the M_D06 groups. 

The median vCG has a mass of 1.7 x 10 13 M Q for both 
M_D06 and DJ306, but with considerable variation. This 
means that the distributions closest to this value are more 
important to the analysis of vCGs than further away. In Fig. 
[9(i0l and [9(b)1 we see that the models are reasonably similar 
at this point, although the peak in the central luminosity 
function is greater in the Munich variant. This suggests that 
the significant lack of vCGs in MJD06 compared to DJ306 
is due to the positions of galaxies in groups rather than the 
absolute numbers, and we refer the reader back to Figure 3. 



6 MAGNITUDE GAP 

The magnitude "gap" between the first- and second- 
ranked (and, indeed, lower-ranked) galaxies within a group 
can be used as a predictor of group (or halo) age 



|von Benda-Beckmann et al.l [2008), as the central galaxy 
tends to grow unceasingly with time via satellite accre- 
tion/stripping. This mechanism inevitably increases the 
magnitude gap. This process is controlled by feedback in 
the central galaxy (§2) and by infalling galaxies. Taken 
to its extreme, such an effect gives rise to the so-called 
"fossil groups", which are groups with a magnitude gap 



cent galaxy infall onto the group (|D'Onghia et al.l 


2005 


Sommer-Larsen 2006; Dariush et al. 2007; Sales et al.l 


2007 


van den Bosch et al. 


20071; Dfaz-Gimenez & Mamonl 


201C 



ever. see lZibetti et al. 20091 ) 



In Fig. 1101 far right column, we compare the magni- 
tude gap distribution between first and second ranked group 
galaxies in the suite of SAMs employed here. We show the 
distribution for all groups (top panel) as well as for the 
LGs, CGs and vCGs separately as we move downwards. The 
M_D06 model (and to a lesser extent, that of M_B07) shows 
a preferred magnitude gap of ~1 mag between the two most- 
luminous galaxies in the model groups, (particularly in the 
CGs and vCGs), while the DJ306 SAM predicts far more 
equal luminosity first- and second-ranked group galaxies. 
The difference in the distributions between the models is 
quite apparent, with a "turn over" in the two Munich mod- 
els, i.e. both are "flatter" and "broader" than the Durham 
models, for all group types and in all mass bins. Again, a sig- 
nificant difference in how the model's galaxies evolve within 
dense environments has been highlighted by these observ- 
able characteristics. 

Dariush et al. (2007; fig 4a) show a comparable rep- 
resentat ion of the to p -left panel of our Fig. 1101 employ- 
ing the ICroton et all l|2006h SAM as applied to the Mil- 
lennium Simulation (in the i?-band), and for a slightly dif- 
ferent mass range, but effectively similar to what we have 
shown. Dariush et al. po int out that the mag nitude gap dis- 
tribution of LGs in the ICroton et ail (2006) model is sim- 
ilar to the In A=2 theoretical model of Milosavlievic et all 
(2006), where In A is the Coulomb logarithm that controls 
the merger rate. When the lCroton et aL I (l2006f) SAM is com- 
pared with the SDSS C4 catalogue (|Miller et all 120051 ) . as 
is done in fig. 4c of Dariush et al,. the apparent mismatch 
at small first- and second-ranked magnitude differences be- 
tween the Munich SAM, and the data become apparent - i.e., 
the SDSS C4 catalogue shows a magnitude difference distri- 
bution which prefers approximately equal luminosity first- 
and second-ranked galaxies in groups and clusters, more con- 
sistent with the Durham SAM predictions. 

If we compare our result to the h alo occupation distri- 
butions of lvan den Bosch et al.l l|2007h . we find that we have 
far fewer fossils groups in the two high mass bins for all 
three models but more for M_D06 groups in the lowest mass 
bin. A halo occupation distribution is a statistical model 
of the number and luminosity of galaxies occupying a dark 
matter halo of a given mass. As such, it is strongly related 
to the conditional luminosity function previously discussed, 
and serve as a base of comparison for the SAMs. 

In order to make a fair comparison of the model pre- 
dictions with the current observations, we take into ac- 
count the selection effects inherent within the data. Specifi- 
cally, the observatio nal result s (i) h ave a limited dynamic 
range of ~2 mag l|Lin et all I1996T ). driven by signal-to- 
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Figure 10. Conditional (mass-dependent) magnitude gaps between first- and second-ranked galaxies for the SAMs included in this 
study. Solid, dashed and dotted lines are for D_B06, M_D06 and M_B07 groups respectively. The top row is for halo-based groups 
divided by group mass and the far right column is for all FoF-groups regardless of mass, specifically, row 2 is for LGs, row 3 is for 
CGs and row 4 is for vCGs. The subsequent panels show the FoF based groups broken down by halo mass. The vertical line shows the 
cut-off for Fossil groups (Sales et al. 2007). The numbers in each panel give the percentage of groups which are fossil systems, and, in 
brackets, the total number of groups in each mass bin. The first number is for D_B06, the second for M_D06 and the third for M_B07. 



noise constraints applied to the lowest luminosity galaxies 
in the survey, and (ii) discard groups that contain fewer 
than four galaxie s within ~2 mag of the first-ranked galaxy. 
iTaeo et all ll200ct) groups were chosen for the comparison be- 
cause the absolute r-band magnitude data was readily avail- 
able. We have imposed comparable selection effects upon 
the models, and the impact upon the luminosity functions 
of the first- and second-ranked CG galaxies is shown in the 
right hand panels of Fig. [11] The left panels show that the 
turnover in the Munich models is no longer apparent, once a 
dynamic range of two magnitudes is imposed upon the (the- 
oretically, infinite) magnitude gap between first- and fourth- 
ranked group galaxies. With these cuts, the models and data 
now lie closer to one another. However, the models produce 
a significant shortage of pairs with low magnitude gaps for 
LGs and a higher population of groups with a magnitude 



gap of one. This effect is more extreme in the Munich mod- 
els but is still present in DJ306. 

What is perhaps more interesting is that the dynamic 
range need only be increased to three magnitudes for the 
models to diverge significantly, with the M_D06 model both 
"broadening" and shifting to lower luminosity, relative to 
the distributions based upon the M_B07 and D_B06 SAMs, 
(right-most column of Fig. I10p . Table 4 shows the popula- 
tions of groups in the two extreme cases of small and large 
magnitude gap for a dynamic range of 3 mag. The differ- 
ences are very large between the Durham and Munich mod- 
els. Certainly, observations with a higher dynamical range 
will provide a good test for differentiating the success of 
the SAMs within group environments. The relative success 
of the different manner that the SAMs implement physical 
processes such as AGN and SN feedback, and how these be- 
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Table 4. Percentage of groups with magnitude gap between first 
and second ranked galaxies greater than 2 magnitudes (top) and 
less than 0.5 magnitudes (bottom) 





gap 


D.B06 


M.D06 


M.B07 


LG 


> 2 


10.8 


20.6 


15.2 


CG 


> 2 


14.9 


28.4 


23.8 


vCG 


> 2 


18.2 


27.5 


29.8 


LG 


< 0.5 


35.6 


22.7 


28.7 


CG 


< 0.5 


30.4 


16.4 


19.7 


vCG 


< 0.5 


27.5 


17.8 


18.7 



come important within group environments where satellite 
accretion modelling is also crucial, can then be better deter- 
mined. One may also ask whether the assumption of sepa- 
rating Central and Satellite populations, whereby only the 
central galaxies experience mergers and no satellite galaxies 
grow while in the group environment, is appropriate when 
two galaxies of almost equal mass often exist within such 
environments. 

The proportion of first-ranked (by luminosity) galaxies 
being centrals is sufficiently high to make the transition from 
the theoretical definitions of "central" and "satellite" galax- 
ies into the observational regime of "brightest" and "second 
brightest" group galaxies - i.e., we can associate the bright- 
est group galaxy with a central, and the second-brightest 
galaxy to a satellite. This then allows us to plot the lumi- 
nosity function of first- (Ml) and second-ranked (M2) group 
galaxies, as shown in the right hand panels of Fig llll and 
associate the distribution in Ml with model centrals, and 
the distribution in M2 with model (brightest) satellites. The 
first-ranked galaxy luminosity function of D J306 is broader 
and flatter than those of the two Munich SAM variants; as 
expected, the M_B07 model galaxies are, on average more 
luminous. For the distribution of second-ranked galaxies, the 
M_D06 galaxies are on average ~1 mag less luminous than 
the Durham model galaxies, and the distribution is broader. 
The right panels highlight that, although the global luminos- 
ity function of galaxies is well matched by observations, the 
distributions for the first and second ranked galaxies shown 
in Fig. 1111 tend to be dimmer and wider than observations. 

In Figure 1121 we demonstrate the impact of imposing 
a dynamic range of three magnitudes between first- and 
fourth-ranked group galaxies; having done so, we find that 
the observations of SDSS4_Y08 match the model predictions 
of DJ306 remarkably well for Loose Groups. This suggests 
that the Durham model, in this regime, provides a better 
match to empirical data than that of the Munich models. 



7 CONCLUSIONS 

By constructing luminosity functions of galaxy groups (rang- 
ing from loose to very compact) using variants of several 
leading SAMs, as applied to the Millennium Simulation, we 
have explored an astrophysical regime in which the SAMs 
have not been inter-compared in gre at detail. Several ob - 
vious differe nces between the M _D06 (|De Lucia et alj|2006l ) 
and DJ306 ([Bower et al.l l200f3 ) (i.e., loosely speaking, the 
Munich and Durham variants, respectively), became ap- 
parent, including an intermediate luminosity "wiggle" in 



the M_D06 group luminosity functions not readily apparent 
when using the D_B06 SAM. We trace the origin of this wig- 
gle to two competing effects resulting from the underlying 
physics within the M_D06 SAM - a steeper faint-end slope 
to the satellite luminosity function and a narrower distribu- 
tion to the central galaxies luminosity function, most likely 
due to the lack of mass stripping in satellite galaxies with- 
out enveloping subhaloes, type 2 groups, and the particular 
formulation of AGN in the Munich models. A systematic 
exploration of parameter space in the respective SAM may 
however, be required to further isolate the cause of the dif- 
ference. 

Observations suggest that such a wiggle in the gro up lu- 
minosity function might exist l|Weinmann et aLll2006l ). sim- 
ilar to that seen when applying the M_D06 SAM. However, 
these same observations tend to show a steeper magnitude 
gap (between first- and second-ranked group members) dis- 
tribution profile than that seen with any of SAMs, and we 
also see significant "flattening" in the M_D06 gap distri- 
bution (i.e., a comparable likelihood for first- and second- 
ranked galaxies to be of equal luminosity, as to have a one 
magnitude luminosity difference), a feat ure that is n o t con - 
sis tent with the data s ets described by iMiller et all (2005) 
or lDariush et all (|2007h . 

The models applied to the Millennium simulation pro- 
duce noticeably different galaxy group properties. The group 
luminosity functions diverge with increasing galaxy density 
meaning that, for example, the cores of clusters in the var- 
ious models have different properties, while the properties 
of the entire cluster will be more similar. As the same dark 
matter background was used in the three models there are 
similar numbers of groups and clusters in the models, but 
according to our definitions, the denser structures are sev- 
eral times more common in the D_B06 model. The MJD06 
model luminosity function shows a peak for the brightest 
galaxies that does not appear in the Durham models and is 
less evident in M_B07. The magnitude gap distributions of 
the models also differ with the Munich and Durham mod- 
els demonstrating a different distribution at the small gap 
part of the distribution. All models show a shallower, wider 
magnitude gap distribution than the observations. This sug- 
gests that improvement in how the central / bright satellite 
luminosities are calculated is required. The designation of a 
single central galaxy which is modelled different manner to 
the other group members is a simplification which may need 
to be improved upon. 

The existence of more, denser, CG and vCG groups in 
the Durham models compared to the Munich sample sug- 
gests that the different merging time-scales and implemen- 
tations of satellite accretion can have noticeable effects on 
the predictions of the models. Similarly the fact that the 
Durham models show a shorter mean galaxy-galaxy sep- 
aration, indicates that these groups are denser. This sug- 
gests that merging time-scales are longer in the Durham 
groups. This is backed up by the luminosity function of 
groups because the evident 'wiggle' in the M_D06 groups 
appears to be due to a smaller population of satellites and 
brighter centrals, which is a direct result of the rate at which 
satellite galaxies are accreted onto the central galaxy. How- 
ever, while observations show a similar 'wiggle' in group lu- 
minosity functi ons, suggesting the s horte r merging time is 
more physical, iMcConnachie et alj l|2008l ) find fewer com- 
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Figure 11. Left column: Distribution of the magnitude gap between 1st and 2nd ranked galaxies, normalised by total number of groups 
in each sample. The thick solid, dashed and dotted lines are for D_B06, M_D06, M_B07 models. The coloured lines show observational 
data. The top left panel shows Tago et al. (2008) groups (red dashed line), SDSS4.Y08 (red solid line) and Tucker et al. 2000 LGs 
(blue dotted line), the middle panel shows the results for CGs and and the bottom left panel shows Allam et al. (2000) compact groups 
(blue dotted line) along with vCGs. Poisson counting uncertainties for each bin are reflected by the accompanying vertical bars. Right 
column: The upper panel shows the luminosity function for second-ranked LG galaxies and the lower panel shows the distribution of 
first-ranked galaxies, and the observational data is taken from Tago et al. (2008) groups (blue line) SDSS4_Y08 (red line). Data here has 
been restricted to mimic a survey in which the magnitude gap between first- and fourth-ranked group galaxies is two. 



pact groups in their field than in the SAMs, suggesting the 
merging time-scale should be even shorter. Contrastingly 
the limited magnitude gap distribution indicates that the 
gap between central and satellite galaxies should be smaller, 
which may be due to additional physics that is not yet im- 
plemented in the models. Our analysis of the timescales of 
merging shows that this is not the satellites in both 

models last a similar amount of time. We emphasises, how- 
ever, that there are more galaxies near the centre of a given 
group/cluster in the Durham models, despite this. We sug- 
gest that this may be due to the additional time the M_D06 
model identifies subhaloes, and reassigns the galaxy posi- 
tion of the central galaxy as the most bound particle in the 
(sub)halo for longer. This may have the effect of keeping 
the galaxy out of the central region for longer. Although we 
do not find a noticeable difference in the merging times of 
galaxies in the two models there is a substantial population 
of groups which do not merge. We can see this because more 
haloes merge in the Durham model but more galaxies merge 
in the Munich models. This serves to build up the number 
of satellites in the cluster, which, fall into the cluster core, 
thus accounting for the observed difference in vCG popula- 
tion and galaxy density distribution. This can explain the 
difference in the magnitude gap distribution because more 



galaxies merge with the central in the Munich model, reduc- 
ing the number of satellites and making the central galaxy 
brighter. 
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Figure 12. Left column: Distribution of the magnitude gap between 1st and 2nd ranked galaxies, normalised by total number of groups 
in each sample. The thick solid, dashed and dotted lines are for D_B06, M_D06, M_B07 models. The coloured lines show observational 
data. The top left panel shows Tago et al. (2008) groups (red dashed line), SDSS4.Y08 (red solid line) and Tucker et al. 2000 LGs 
(blue dotted line), the middle panel shows the results for CGs and and the bottom left panel shows Allam et al. (2000) compact groups 
(blue dotted line) along with vCGs. Poisson counting uncertainties for each bin are reflected by the accompanying vertical bars. Right 
column: The upper panel shows the luminosity function for second-ranked LG galaxies and the lower panel shows the distribution of 
first-ranked galaxies, and the observational data is taken from Tago et al. (2008) groups (blue line) SDSS4_Y08 (red line). In contrast 
with Fig. 1111 data here has been restricted to mimic a survey in which the magnitude gap between first- and fourth-ranked group galaxies 
is three. 
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